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ABSTRACT 

Recently a hard X-ray transient event, Sw J1644+57, was discovered by the Swift satellite. It likely 
marks the onset of a relativistic jet from a supermassive black hole, possibly triggered by a tidal 
disruption event. Another candidate in the same category, Sw J2058+05, was also reported. The low 
event rate suggests that only a small fraction of TDEs launch relativistic jets. A common speculation 
is that these rare events are related to rapidly spinning black holes. We attribute jet launching to the 
Blandford-Znajek mechanism, and use the available data to constrain the black hole spin parameter 
for the two events. It is found that the two black holes indeed carry a moderate to high spin, suggesting 
that black hole spin is likely the crucial factor behind the Sw J1644+57 - like events. 
Subject headings: accretion, accretion disks-black hole physics-magnetic fields 



1. INTRODUCTION 

A hard X-ray transient event, Swift J16449. 3+573451 
("Sw J1644+57" hereafter) was discovered by the Swift 
satellite (Gehrels et al. 2004) on 2011 March 28, initially 
reported as GRB 110328A (Cummings et al. 2011). It 
was soon realized that it is not a regular GRB. Long-term 
follow up observations with Swift (Burrows et al. 2011) 
revealed that it has extended emission in the X-ray band 
without significant decay. A much longer variability time 
scale St ~ 100 s than normal GRBs (Burrows et al. 2011) 
as well as its location near the center of a z — 0.354 host 
galaxy (Levan et al. 2011) link the source to a super- 
massive black hole (BH) with M, ~ 1O 7 M (Burrows 
et al. 2011; Bloom et al. 2011). The super-Eddington 
X-ray luminosity (Burrows et al. 2011), bright radio af- 
terglow (Zauderer et al. 2011), as well as a historical 
stringent X-ray flux upper limit suggest that the event 
marks the onset of a relativistic jet from a super-massive 
BH. The sharp onset, and gradual fade-away of X-ray 
flux suggest that the transient may be triggered by tidal 
disruption of a star by the BH (Bloom et al. 2011; Bur- 
rows et al. 2011). 

A second candidate in the same category, Swift 
J2058.4+0516 ("Sw J2058+05" hereafter) was reported 
by Cenko et al. (2011), who argued that its observational 
properties are rather similar to those of Sw J1644+57. 

These events are rare (2 detected in 7 years of Swift 
operation). Most previous studies of tidal disruption 
events (TDEs) from super-massive BHs did not expect 
an associated relativistic jet (e.g. Kobayashi et al. 2004; 
Strubble & Quataert 2010; Lodato & Rossi 2011, but 
see Lu et al. 2008; Giannios & Metzger 2011). The 
neutrino-annihilation mechanism as invoked in GRBs is 
inadequate given a much lower accretion rate as com- 
pared with GRBs (e.g. Shao et al. 2011). A plausible 
mechanism to launch a jet from such an event is to tap 
spin energy of the BH through a magnetic field, which 
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connects the BH event horizon and a remote astrophysi- 
cal load (Blandford & Znajek 1977, here after BZ). This 
mechanism has been widely invoked to interpret AGN 
jets (e.g. Vlahakis & Konigl 2004). Within this scenario, 
a dominant magnetic jet composition is envisaged. In- 
deed, modeling the emission of Sw J1644+57 suggests 
that the jet is highly "particle starved" , i.e. most energy 
is not carried by particles. A natural inference is that 
the outflow is Poynting-flux-dominated (Burrows et al. 
2011). 

If one identifies the BZ mechanism as the jet launch- 
ing mechanism, the BH spin parameter (a dimensionless 
angular momentum of the BH, which ranges from for 
no spin to 1 for the maximum spin) can be constrained 
from the data. This is the subject in this Letter. 

2. METHOD 

The bipolar BZ jet power from a BH with mass M, 
and angular momentum J is (Lee et al. 2000; Li 2000; 
Wang et al. 2002; McKinney 2005) 



L BZ = 1.7 x 10 
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where a. = Jc/(GM%) is the BH spin parameter, M #) i 
M./1O 6 M , B. fi = S./10 6 G and 

F(a.) = [(l + q 2 )/q 2 }[(q + l/q) arctanq - 1] 



here q = a./(l + ^1 - of), and 2/3 < F{a,) < tt - 2 for 
< a. < 1. It apparently depends on M # , B t , and a,. 

However, since an isolated BH does not carry a magnetic 
field, B, is closely related to the accretion rate M and 
the radius of the BH, which depends on M,. Combin- 
ing these dependences, one finds that Lbz is essentially 
independent of M,, but is rather a function of M and 
a.. This may be proven with the following rough scal- 
ings. For a Newtonian disk, angular momentum equa- 
tion states Mr 1 {GM,lr' i ) 1 l 2 ~ — 47rr 2 r r( ^/i, where T r( j, 
is the viscous shear, and ft oc r is the half disk thick- 
ness, with h/r -C 1 for a thin disk and h/r ~ 1 for a 
thick disk. Adpoting the a-prescription for viscosity, the 
viscous shear can be expressed as T r( f> — ~aP, where 
P is the total pressure in the disk. One therefore de- 
rives P ~ Mr 2 (G?M./r 3 ) 1 / 2 /(47rr 2 a/i) cx (M./r 5 ) 1 / 2 cx 



2 



Af. 2 , where we have applied the scaling r oc r s oc Af. 
(r s = 2GM,/c 2 is the Schwarzschild radius). To make 
an efficient BZ jet, the accretion inflow should carry a 
large magnetic flux (e.g. Tchekhovskoy et al. 2011). It 
is reasonable to assume that magnetic fields in the disk 
are in close equilibrium with the total pressure, so that 
P 2 oc P oc Af~ 2 . Inserting this dependence to Eq.([T]), 
one finds Lbz is essentially independent of Af.. 

More precisely, we adopt the following prescription to 
treat the problem. The total pressure in the disk can be 
expressed as P = P ra( j + P gas + Pb, where P ra( j = aT 4 /3, 
Pgas = nkT, and Pb = B^ isk /8ir are radiation, gas and 
magnetic pressure, respectively. Here T is the tempera- 
ture of the disk, n is the gas particle number density, a 
is radiation density constant, and k is Boltzmann con- 
stant. We denote Pb = /3P, and take j3 ~ 0.5 in this 
work. This corresponds to a maximized magnetic flux. 
A smaller ft would demand an even larger a. to reach 
a same BZ power. So taking f3 ~ 0.5 gives an estimate 
of a conservative lower limit of a.. For the two sources 
(Sw J1644+57 and Sw J2058+05) we are interested in, 
the accretion rate is estimated close to or even larger 
than the Eddington accretion rate (Burrows et al. 2011; 
Cenko et al. 2011). In this regime, a thin disk model 
may be adequate to describe the disk, and we apply it to 
estimate disk pressure for simplicity. The disk pressure 
peaks in the inner region where radiation pressure may 
dominate. As a result, for (3 ~ 0.5, one has 
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where the temperature of the disk (a.- and r-dependent) 
can be written 
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where a is the Stephan-Boltzmann constant, and the gen- 
eral relativistic correction factor (Page & Thorne 1974) 
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where \ 
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(fin/r 9 ) 1/2 , and r g 



GAf./c 2 . For a Kerr BH, the disk inner edge r ln is ex- 
pressed as (Bardeen et al. 1972), 



r m /r g = 3 + Z 2 - [(3 - Z x ){2> + Z± + 2Z 2 )] 
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for < a. < 1, where Zi = 1 + (1 - a 2 ) 1 / 3 [(l -I- a.) 1 / 3 + 
(1 - a.) 1/3 ], Z 2 = (3a 2 + Ziyi*. In Eq.©, Xi,X2,X3 
are the three roots of X 3 — 3x + 2a. = 0, i.e. xi = 
2 cos(| cos -1 a. — 7r/3), X2 = 2 cos(| cos -1 a. + 7r/3), 
X.3 = — 2cos(icos~ 1 a.). It is easy to check that f(r = 



ri n ) = and /(r» ?'j n ) ~ 1. For a Newonian disk, / can 
be simply written as / = 1 — y/r^Jr. 

We assume that the strength of the magnetic field 
threading the BH is comparable to the largest field 
strength in the disk, i.e. B, ~ PJfJsk- For given M and 
a,, we solve numerically the temperature profile of the 
thin disk and identify the radius r pca k where T reaches 
the maximum. It is found that r poa k is very close to 
ri n , the inner most radius of the accretion disk. We then 
calculate from Eq.([3]), which is assigned to P.. Ap- 
plying Eqs.([T|) and ([2]), one then obtains Lbz once M and 
a, are specified. 

Observationally, a time-dependent isotropic X-ray lu- 
minosity Lx.iso w as measured for the two sources. This 
luminosity can be connected to the BZ power through 



vLbz = fbLx,: 



(7) 



where r\ is the efficiency of converting BZ power to X-ray 
radiation, 



fb = 



An 

An 




< 1 



(8) 



is the beaming factor of the jet (Burrows et al. 2011), 
Af2 is the solid angle of the bipolar jet, and T and 9j are 
the Lorentz factor and opening angle of the jet. 

In this work, we adopt 77 ~ 0.5, motivated by a 
potentially high radiation efficiency of a magnetically- 
dominated jet (e.g. Drenkhahn & Spruit 2002; Zhang & 
Yan 2011). Again this gives a conservative lower limit of 
the inferred a,, since a less efficient jet would demands 
an even higher spin rate in order to interpret the same 
observed luminosity. 

The parameter fb has a large uncertainty. Based on 
the current data, one cannot precisely measure T and 
9j. Here we apply a statistical method to estimate 
the range of fb within the TDE framework (see also 
Burrows et al. 2011). First, the facts that two such 
events (Sw J1644+57 and Sw J2058+05) were detected 
by Swift in ~ 7 years and that the field view of Swift 
Burst Alert Telescope (BAT, Barthelmy et al. 2005) is 
~ 47r/7 sr suggest that the all-sky rate of such events 
is P bs ~ 2yr _1 , with a 90% confidence interval of 
(0.44 - 5.48)yr~ 1 (Kraft et al. 1991). Next, the TDE 
rate is estimated as ~ 10 -5 — 10~ 4 yr~ 1 galaxy^ 1 based 
on observational (Donley et al. 2002; Gezari et al. 2009) 
theoretical (Wang & Merritt 2004) studies. The galaxy 
number density is n ga i ~ 10~ 3 — 10~ 2 Mpc~ 3 (Tundo et 
al. 2007). Sw J2058+05 was marginally detected at 
z = 1.1853 (Cenko et al. 2011). We then obtain the 
total TDE event rate within the volume (z < 1.1853) 
i? tot - 10 4 - 10 5 yr _1 - Considering that only ~ 10% of 
the population can launch a jet (the "radio-loud" AGN 
fraction, Kellerman et al.1989; Cirasuolo et al. 2003), 
the beaming factor can be estimated as 



fb 



Robs 



10%P t 



G (4.4 x 10~ b ,5.5 x 10~ 3 ) 



(9) 



Finally, in order to infer a. from Lx.isoi fb (given 
ft = 0.5 and r\ = 0.5), one needs to know the accretion 
rate Af. This is an even loosely constrained parameter. 
However, if one assumes that the luminosity history of 
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the light curve well delineates the accretion history of the 
BH (noticing Lbz cx M), one can normalize the accretion 
rate using the total accreted mass based on the observed 
flux and fluence of the event. For Sw J1644+57, since 
the source has entered a decaying phase, during which 
the residual fluence no longer contributes significantly to 
the total fluence, one can take the current total X-ray 
fluence as a good proxy of the total mass of the tidally 
disrupted star. Taking the peak flux as example, the 
peak accretion rate can be estimated as 



M pea k = 



i^ cak (l + z) 



M* = 



L 



peak 



E- 



X,iso 



(10) 



where i7 , P eak j s the peak X-ray flux, Sx is the total X- 
ray fluence, 1^ is the peak isotropic X-ray luminosity, 
Ex.iso is the isotropic X-ray energy, and M* is the to- 
tal mass of the star. The factor (1 + z) was applied to 
convert the observed time to the time in the source rest 
frame. The accretion rate at other epochs can be esti- 
mated similarly. The range of M* may be between O.IMq 
and IOMq. One can then derive a mass-dependent con- 
straint on a.. 

3. SW J1644+57 AND SW J2058+05 

Now we apply the above method to the two sources. 

According to Burrows et al. (2011), during the first 
50 days after the first BAT trigger the total X-ray en- 
ergy corrected for live-time fraction for Sw J1644+57 is 
-Ex,iso(J1644) - 2 x 10 53 erg in the 1.35-13.5 keV rest- 
frame energy band (corresponding the energy band of 
Swift X-Ray Telescope, Burrows et al. 2005). The peak 
luminosity in the same energy band is L x ca 3 k 3 (J1644) ~ 

2.9 x 10 48 erg s _1 . The accretion rate (Eq.lfTD))) can be 
estimated as 

M poak (J1644) ~ 1.45 x 10- 5 M*s _1 . (11) 

Fig UK presents the constraint on a, for Sw J1644+57 
as a function of M„, with /3 = 0.5, r\ — 0.5. The two 
boundary lines correspond to two ends of the range of 
/b, with the lower and upper boundaries correspond to 
fb = 4.4 x 10~ 5 and f b = 5.5 x 10 -3 , respectively. The 
middle dashed line corresponds to the most probable 
value ft, ~ 10~ 3 . One can see that the supermassive BH 
is demanded to have a moderate to high spin rate. Given 
the standard stellar initial mass function, the number of 
low-mass stars is much more abundant than the high- 
mass stars. If one takes M* = 1M Q , the required range of 
a. is (0.23, 0.85), with the most probably value a. = 0.63. 
For M* = O.IMq (more probable), the range of a, is 
(0.51, 0.98) with the most probable value a, = 0.90. 

Sw J2058+05 was discovered by Swift BAT through 
a 4-day (2011 May 17-20) integration. A subsequent 
target-of-opportunity (ToO) observation 8 days after the 
end of 4-day integration (2011 May 28) revealed an X- 
ray source that behaves very similarly to Sw J1644+57 
(Cenko et al. 2011), suggesting that it is very likely an- 
other Sw J1644+57-like event. The 0.3 - 10 keV peak 
flux is Fxlso — x 10~ n erg cm' 2 s" 1 , corresponding 
to a peak luminosity L^f ~ 3 x 10 47 erg s _1 . The total 
X-ray fluence from the beginning of the ToO observation 
to 2011 July 20 is S x ^ 1.0 x 10~ 4 erg cm" 2 (Cenko et 



al. 2011). Since we did not catch the X-ray emission 
during the first 11 days (May 17 - 27) when it is sup- 
posed to be much brighter, the registered X-ray fluence 
only corresponds to a small fraction of the total mass of 
the star. One can still apply Eq.(fT0|) to estimate the ac- 
cretion rate, except that one should replace M* by CM,, 
where £ < 1 is the fraction of stellar mass that is accreted 
after May 28. Noticing z = 1.1853 (Cenko et al. 2011), 
this gives a peak accretion rate 



M P cak(J12058) ~ 1.72 x Kr 7 C-iM»s -1 



(12) 



where £ = 0.1£_i has been adopted. 

The constraint on BH spin for Sw J2058+05 is shown 
in FigQ}). We find that the demand for BH spin is even 
more stringent for this source. For M* = 1M Q , the re- 
quired range of a, is (0.49, 0.98), with the most probably 
value a, = 0.89. For M* = O.IMq (more probable), the 
range of a, is (0.81, 0.998) with the most probable value 
a. = 0.99. 

4. CONCLUSION AND DISCUSSION 

Sw J1644+57 and Sw J2058+05 are the first two proto- 
type objects in the newly identified astrophysical phe- 
nomenon, namely, a relativistic jet associated with a 
TDE from a supermassive BH. A straightforward ques- 
tion is why only some TDEs launch jets. Based on ob- 
servational properties, it has been argued that the jet 
has to be Poynting-flux dominated (Burrows et al. 2011; 
Shao et al. 2011). Invoking the BZ mechanism as the 
power of the jet, we show here that both events need 
to invoke a BH with a moderate to rapid spin in order 
to interpret the observations: the most probably values 
are <z.(J1644) = 0.63,0.90 and a.(J2058) = 0.89,0.99 for 
M* = 1, O.IMq, respectively. We therefore suggest that 
BH spin is the key factor behind the Sw J1944+57 - like 
events, although other factors may also play a role (e.g. 
Cannizzo et al. 2011). 

An elegant feature of the method we employ is that 
the inferred BH spin parameter a. essentially does not 
depend on the BH mass. On the other hand, knowing a. 
leads to better constraint on BH mass based on the vari- 
ability argument. For example, the observed minimum 
variability time scale of J1644+57 is <5£ bs,min ~ 100 s. If 
one relates <5t m ; n = #t bs,mm/ (1 + z ) to the time scale de- 
fined by the innermost radius of the accretion disk, ri n /c, 
one can derive the BH mass of Sw J1644+57 



AT.., 



15| ^ 

r„ 



St 



obs 



100s 



(13) 



One has 2.5 < Mg < 15 for < a, < 1, with the most 
probably value M§ = 6.5 for a. = 0.9. This is consistent 
with the constaincd BH mass from the M — Lbuige rela- 
tion, which gives an upper limit of 2 x 1O 7 M (Burrows 
et al. 2011). 

From Eq.fl]), one can also infer the strength of the 
magnetic field at the BH horizon 



B,, 



(14) 



Taking BH mass M. = 6.5 x 1O 6 M and the most probaly 
value for BH spin a. = 0.9, one finds that the magnetic 
field threading BH would be B. ~ 1.1 x 10 6 G, which 
is much higher than the average field strength of a typ- 
ical main sequence star (< 10 3 G). The accumulation 
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of magnetic flux by accretion flow and instability in the 
disk may account for such high magnetic field stength 
(e.g. Tchekhovskoy et al. 2011). 

For Sw J2058+05, due to the low X-ray flux at the late 
epochs a much looser constraint on variability (5i bs.min < 
10 4 s was obtained (Cenko et al. 2011), so that the pre- 
cise values of M, and B, cannot be derived. 

An alternative scenario to interpret Sw J1644+57 - like 
event may be the onset of an AGN (Burrows et al. 2011). 
This scenario, which predicts that the two sources will 
be active at least in the following millennium, may be 
less favored in view of the rapid onset of emission in Sw 
J1644+57 and the gradual decay in both events, but is 
not ruled out. Our method can be applied to this sce- 
nario as well, except that fb becomes much larger (due 
to the intrinsic rarity of AGN onset events), and Af* is 
no longer limited to the range of 0.1 — 10M Q . Our anal- 
ysis suggests that the BZ power is likely not adequate to 
interpret the data (because of the large emission power 
demanded by the small /& factor) even for maximum spin 
(a. ~ 1), unless the accretion rate is much higher, so that 
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the total amount of fuel M* 3> IMq. This is not impos- 
sible since the fuel on the AGN onset scenario is from a 
gas cloud near the BH, whose mass is not specified. 

Krolik & Piran (2011) proposed a model for Sw 
J1644+57 invoking a white dwarf being tidally disrupted 
by a smaller BH (Al, ~ 10 4 M Q ). Regardless of how 
this model may interpret the <5i m i n and the apparent as- 
sociation of the source with the center of host galaxy, 
the derived a, range also applies to their model for the 
M* range of a white dwarf, since our constraint is M.- 
independent. 

Finally, in our calculations we did not consider evolu- 
tion of a. during the accretion phase. This is justified 
given the large M./M* ratio. 
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Fig. 1. — Parameter space of M* — a. for Sw J1644+57 (left panel) and Sw J2058+05 (right panel) with /3 = 0.5 and r\ = 0.5. The shaded 
region indicate the range for BH spin bracketed by ft, = 4.4 X 10~ 5 (lower) and ft = 5.5 X 10~ 3 (upper). The dashed line corresponds to the 
most probable value fi = 10 . For Sw J2058+05, C, = 0.1 is adopted. It is shown that the most probable values are a, (J1644) = 0.63, 0.90 
and a.(J2058) = 0.89,0.99 for M* = 1,0.1M Q , respectively. 



